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ABSTRACT 

We report on the detection of X-ray emission from the unusual Galactic radio source G357. 7-0.1 (the "Tor- 
nado"). Observations made with the Chandra X-ray Observatory demonstrate the presence of up to three sources 
of X-ray emission from the Tornado: a relatively bright region of dimensions 2' x 1' coincident with and interior 
to the brightest radio emission at the "head" of the Tornado, plus two fainter extended regions possibly associated 
with the Tornado's "tail". No X-ray point sources associated with the Tornado are seen down to a 3er luminosity 
(0.5-10 keV) of 1 x 10 33 erg s" 1 , for a distance to the system of 12 kpc. The spectrum of the brightest region of 
X-rays is consistent with a heavily absorbed (Nh ~ 10 23 cm" 2 ) thermal plasma of temperature kT ~ 0.6 keV; an 
absorbed power law can also fit the data, but implies an extremely steep photon index. From these data we tenta- 
tively conclude that the Tornado is a supernova remnant (SNR), although we are unable to rule out the possibility 
that the Tornado is powered either by outflows from an X-ray binary or by the relativistic wind of an unseen pulsar. 
Within the SNR interpretation, the head of the Tornado is a limb-brightened radio shell containing centrally-filled 
thermal X-rays and which is interacting with a molecular cloud. We therefore propose that the Tornado is a 
"mixed morphology" supernova remnant. The unusual tail component of the Tornado remains unexplained in this 
interpretation, but might result from expansion of the SNR into an elongated progenitor wind bubble. 

Subject headings: ISM: individual (G357. 7-0.1) — ISM: supernova remnants 



1. INTRODUCTION 

The bright radio source G357. 7-0.1 is extended, has high 
levels of polarization and shows a non-thermal spectrum, prop- 
erties which at first led to its classification as a supernova 
remnant (SNR). However, high-resolution images revealed a 
bizarre axial "head-tail" structure (Shaver et al. 1985; Helfand 
& Becker 1985; top panel of Fig 1), for which this source was 
nicknamed "the Tornado". The morphology of the Tornado 
is very different from the approximately circular morphologies 
seen in other SNRs. This has led to various claims that the Tor- 
nado is an exotic SNR expanding into an unusual environment, 
is a nebula powered by a high velocity pulsar, or results from 
precessing jets from an X-ray binary (e.g. Manchester 1987; 
Shull, Fesen, & Saken 1989; Stewart et al. 1994). 

Recent observations have provided some new clues. Frail 
et al. (1996) detected a 1720-MHz OH maser from the "head" 
of the Tornado. The 1720 MHz emission, when seen with- 
out the other OH maser lines, is produced by shocks propagat- 
ing into molecular clouds; further evidence for an interaction 
with molecular material is provided by the presence of emission 
from 13 CO 1-0 and shock-excited 2.12-/xm H2 at this same po- 
sition and velocity (Lazendic et al. 2003). The maser velocity, 
also now confirmed by H I absorption measurements (Brogan 
& Goss 2003), implies a distance to the Tornado of 12 kpc. 

Despite these new data, the Tornado still defies classifica- 
tion. Most problematic has been the failure to detect emis- 
sion in other wavebands. We consequently examined archival 
ASCA data of this region, and identified a weak (~ 3a) X-ray 
source coincident with the Tornado's head; this marginal detec- 



tion was also briefly noted by Yusef-Zadeh et al. (2003). We 
here present an observation of the Tornado with the Chandra 
X-ray Observatory, aimed at confirming this source of X-rays, 
and studying its morphology and spectrum. 

2. OBSERVATIONS AND RESULTS 

The Tornado was observed with Chandra on 2002 May 07. 
We used the ACIS-I detector, consisting of an array of four 
front-illuminated CCDs and giving a field of view approxi- 
mately 17' x 17'. The aimpoint, located near the center of the 
array, corresponded to coordinates RA (J2000) 17 h 40 m 15?33, 
Decl. (J2000) -30°57'59."2. The effective exposure time was 
20723 seconds, and was free of background flares. 

2.1. Imaging 

Data were analyzed using CIAO v2.3. We used the level-1 
events supplied by the Chandra X-ray Center (CXC), and then 
corrected them for charge-transfer inefficiency using standard 
CXC routines. An image was then formed using events in the 
energy range 0.3 to 10 keV; point sources in this image were 
identified using the wavdetect algorithm within CIAO. To 
highlight diffuse emission, we applied the adaptive smoothing 
algorithm csmooth. We corrected for off-axis vignetting and 
for the variable exposure produced by chip gaps and bad pixels 
by forming an exposure map, which we then smoothed with 
the same kernel as was applied to the image. Applying the 
smoothed exposure map to the smoothed image resulted in a 
flux-calibrated image of diffuse emission. The resulting X-ray 



image is shown in the lower panel of Figure 1 , demonstrating 
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X-RAY EMISSION FROM THE TORNADO 



that faint emission is detected from three distinct regions coin- 
cident with the Tornado, which we label A, B and C. 

The brightest X-rays are from region A, lying near the "head" 
of the Tornado, and comprising a ~ 2' x 1' elliptical region 
elongated along a position angle ~ 150° (north through east). 
Comparison with radio data shows that these X-rays are wholly 
contained within the head region, occupying the western half of 
the head's interior. While region A is brighter and more con- 
centrated at its south-eastern tip, no unresolved source is iden- 
tified in this region by wavdetect, nor is apparent by visual 
inspection of the unsmoothed image. 

Regions B and C represent two other possible sources of 
extended X-rays from the Tornado, located <~ 2' east of the 
head region. Region B is a very faint circular clump of ex- 
tent ~ 30", while region C is a brighter and more elongated 
feature of extent ~ 1'. Given the faintness of regions B and 
C, we simply note their existence, and defer detailed studies of 
these regions until deeper observations have been performed. 
The wavdetect algorithm identifies two X-ray point sources 
lying within the radio boundaries of the Tornado, as marked in 
Fig 1 . Both sources are of marginal significance, each contain- 
ing only 5-6 counts. 

2.2. Spectroscopy 

Of the three extended regions of X-ray emission identified 
in §2.1, only region A is bright enough from which to obtain 
a useful spectrum. Events were extracted from an elliptical 
extraction region of dimensions 150" x 120" centered on this 
emission, and the data (source plus background) were then re- 
grouped so that there were at least 50 counts in each energy 
bin. The effective area for this source spectrum was determined 
by weighting by the brightness of the source over the extrac- 
tion region using the CI AO task mkwarf. The local back- 
ground was determined by extracting a spectrum from an adja- 
cent emission-free region of the CCD immediately to the south- 
west of region A. 

The corresponding spectrum contains 162 ± 37 background- 
subtracted counts, which fall into 12 bins spread over the energy 
range 0.3-10 keV. We fit two simple models to this spectrum: a 
power law, and a Raymond-Smith plasma in collisionally ion- 
ized equilibrium, 6 both modified by foreground photoelectric 
absorption. Both models give good fits to the data, the parame- 
ters for which are listed in Table 1 . However, the implied pho- 
ton index for the best power law fit is unusually steep, Y > 4.5 
at 90% confidence. In Table 1 we also present power law fits 
for fixed photon indices Y = 2.3 and Y = 3.5 (representing ex- 
treme value for observed sources; see §§3.2 & 3.3), which give 
comparatively poor fits to the data. 

3. DISCUSSION 

We have clearly identified X-ray emission from the head of 
the Tornado; the observed count-rate is consistent (to within 
~ 50%) with the level of emission suggested from archival 
ASCA data. Although we detect relatively few X-ray photons 
from the Tornado, these data provide new constraints on the 
nature of this complicated source. 

3.1. Outflows from an X-ray binary? 

Stewart et al. (1994) propose that the Tornado is powered by 
twin outflows from a central X-ray binary, presumably located 
in the center of and along the symmetry axis of the overall radio 

6 More complicated thermal models including non-equilibrium ionization give 



nebula, and which generate opposed precessing jets as is seen 
for SS 433 in the SNR W50 (e.g. Margon 1984). We argue in 
§3.3 below that region A has a thermal spectrum. In the context 
of the situation proposed by Stewart et al. (1994), these thermal 
X-rays might represent the working surface where a collimated 
outflow from the central source collides and interacts with sur- 
rounding gas, as is seen at the termination of the eastern lobe 
of W50 (Brinkmann, Aschenbach, & Kawai 1996; Safi-Harb & 
Ogelman 1997). 

Of immediate concern for such a model is that no bright 
unresolved X-ray emission, corresponding to a central pow- 
ering source, is seen along the axis of the Tornado. The 3a 
upper limit on the count-rate from any such source is « 5 x 
10~ 4 cts s" 1 in the energy range 0.3-10 keV. Assuming a fore- 
ground absorbing column density of Nh ~ 1 x 10 23 cm" 2 (see 
§3.3 below) and a power-law spectrum with photon index Y = 2, 
we infer an unabsorbed flux (0.5-10 keV) for such a source 
/x^6x 10~ 14 ergs cm" 2 s" 1 , which implies an upper limit on 

the X-ray luminosity (0.5-10 keV) L x <J 1 x 10 33 ergs s" 1 at a 
distance of 12 kpc. 

From this non-detection an exact analogy with SS 433 can 
immediately be ruled out: even when not in outburst, such a 
source would have an X-ray luminosity Lx ~ 10 36 ergs s" 1 and 
a 1.4-GHz flux density ~ 100 mJy (see Margon 1984 and refer- 
ences therein), which would have been easily detected in Chan- 
dra and VLA observations of the Tornado, respectively. How- 
ever, we below note two other X-ray binaries known to gener- 
ate outflows, which if embedded in the Tornado might not be 
directly detectable. 

The recently discovered high-mass X-ray binary LS 5039 
generates steady jets of velocity ~ 0.2c and kinetic luminosity 
~ 10 37 ergs s" 1 (Paredes et al. 2000), which is consistent with 
the rate of energy input required to power the radio emission 
seen from the Tornado (Helfand & Becker 1985). The X-ray 
luminosity of LS 5039 possibly varies as a function of orbital 
phase, but at its minimum is only slightly higher than the upper 
limit determined for any point source in the Tornado (Reig et al. 
2003). At a distance of 12 kpc, the 1.4-GHz flux density of this 
source would only be 1-3 mJy (Marti, Paredes, & Ribo 1998), 
too faint to be seen in existing VLA radio images of this region. 

The low-mass X-ray binary 4U 1755-338 is currently in a 
quiescent state with L x < 4 x 10 31 ergs s" 1 , but previously un- 
derwent an active phase lasting several decades, in which it gen- 
erated X-ray jets ^4 pc in extent (Angelini & White 2003). The 
Tornado could similarly have been powered by a period of pro- 
longed activity, but is now quiescent. Jets of the surface bright- 
ness seen for 4U 1755-338 would be too faint to be detected 
here. 

We note that an alternative X-ray binary model has been pro- 
posed by Helfand & Becker (1985), who argue that the Tornado 
is powered by source at the very western edge of the head, and 
that the high velocity of the binary then gives the Tornado its 
unusual appearance. The arguments invoked above can explain 
the failure to detect X-ray emission from the binary itself. How- 
ever, in this model the energetic interaction between the outflow 
and its surroundings should occur only at the very western edge 
of the Tornado. In such a situation, there is no simple explana- 
tion for the presence of extended X-ray emission from region 
A, since this should be neither a region of active particle accel- 
eration nor an area in which shock-heated gas is produced. We 
therefore think this possibility unlikely. 

■ identical best-fit parameters. 



GAENSLER ET AL 



3 



3.2. A pulsar powered nebula? 

Shull et al. (1989) have argued that the Tornado is a pulsar 
wind nebula (PWN) powered by a high-velocity pulsar. In such 
a situation, the pulsar is presumed to be inside of or to the west 
of the head of the Tornado. We can consequently expect to 
see extended synchrotron X-ray emission from a surrounding 
PWN, plus possibly unresolved magnetospheric X-ray emis- 
sion from the pulsar itself. 

Region A clearly corresponds to X-ray emission in the ex- 
pected location for this interpretation, and its spectrum can be 
well-fitted by a power law characteristic of synchrotron emis- 
sion. However, all known PWNe show power law X-ray spec- 
tra with photon indices in the range 1.3 <; T <^ 2.3 (Gotthelf 
2003); Table 1 shows that a photon index even at the upper 
limit of this range, T = 2.3, gives a comparatively poor fit, 
which we can exclude at the ~ 2.3er level. While more data 
are needed to confirm whether this source indeed has an anoma- 
lously steep spectrum, this measurement suggests that the prop- 
erties of region A are not consistent with those expected for a 
pulsar-powered nebula. We note that while we detect no point 
source embedded in region A, the corresponding upper limit 
L x £ 1 x 10 33 ergs s" 1 (see §3.1 above) is above the X-ray lu- 
minosities of many young pulsars (Possenti et al. 2002). 

3.3. A supernova remnant 

The Tornado shows bright, polarized, limb-brightened radio 
emission with a non-thermal spectrum, properties which re- 
sulted in its original classification as a SNR. However, one must 
then explain its bizarre and unique morphology. 

In a dense ambient medium a supernova progenitor's space 
velocity can take it near the edge of, or even outside of, its cir- 
cumstellar bubble before the star explodes. The SNR resulting 
from such a system will have a highly elongated axial morphol- 
ogy (Rozyczka et al. 1993; Brighenti & D'Ercole 1994), re- 
sembling that seen for the Tornado. The OH maser identified 
towards this source by Frail et al. (1996), along with the recent 
detection of shock-excited 2.12-/im H2 emission and surround- 
ing 13 CO emission by Lazendic et al. (2003), all argue that 
indeed the Tornado is expanding into a dense environment. 

For a SNR interpretation there are several different possi- 
bilities for the X-ray emission expected. We might see limb- 
brightened synchrotron emission from shock-accelerated elec- 
trons (e.g. G347. 3-0.5; Slane etal. 1999), limb-brightened ther- 
mal emission from shock-heated ambient gas (e.g. IE 0102.2- 
7219; Hughes, Rakowski, & Decourchelle 2000), or centrally- 
filled thermal emission from hot gas in the interior, correspond- 
ing to a so-called "mixed morphology" or "thermal composite" 
SNR (e.g. W44; Rho et al. 1994). 

The X-ray photon index for synchrotron-emitting SNRs is in 
the range 2.5 <; T <; 3.5 (see Petre 2001 for a review) which, as 
for the case of a PWN in §3.2, is distinctly flatter than the best- 
fit power law observed. On the other hand, the temperature of 
kT ~ 0.6 keV implied by thermal models is quite reasonable for 
thermally emitting SNRs, as is the corresponding unabsorbed 
X-ray luminosity (0.5-10 keV) of L x ~ 10 36 erg s" 1 for a dis- 
tance of 12 kpc. The foreground column Nh ~ 10 23 cm" 2 is 
high, but is well below the total integrated column in this direc- 
tion of N H ~ 1.5 x 10 23 cm" 2 (Dickey & Lockman 1990). 

If we identify the head of the Tornado as a limb-brightened 
radio SNR, then the emission seen with Chandra originates 
wholly from the SNR interior, and we correspondingly can 
tentatively classify the Tornado as a mixed morphology SNR 



(as had been earlier speculated by Yusef-Zadeh et al. 2003 
from the possible presence of ASCA emission). Whether such 
centrally-filled thermal X-rays are produced by evaporation of 
dense clouds overrun by the shock, through thermal conduction 
in the interior, or via some other process is currently controver- 
sial (Rho & Petre 1998). However, regardless of the mechanism 
invoked, there have been repeated arguments that such a SNR is 
the result of an interaction with a dense molecular cloud (Green 
et al. 1997; Rho & Petre 1998; Yusef-Zadeh et al. 2003). The 
OH maser, 2.12-/im H2 and 13 CO emission seen here all sup- 
port this interpretation. 

Wardle (1999) has pointed out that the soft X-ray emis- 
sion produced by a mixed morphology SNR can produce the 
ionizing flux ( ;> 10" 16 s" 1 needed to disassociate H2O and 
produce the OH seen in maser emission. At a distance of 
12 kpc, the maser produced by the Tornado is <~ 3.5 pc from 
the source of X-rays, implying a flux as seen by the OH maser 
fx ~ 6 x 10" 4 ergs cm" 2 s" 1 . The consequent ionizing flux is 
( = N e af x ~3x 10" 15 s" 1 , where N e w 30 keV" 1 is the number 
of electrons produced per keV of ionization and a « 2.6 x 10~ 22 
is the photoabsorption cross-section of hydrogen at 1 keV (War- 
die 1999). This rate of ionization is more than sufficient to con- 
vert water into OH molecules. 

4. CONCLUSIONS 

A short observation of G357. 7-0.1 with the Chandra X- 
ray Observatory has confirmed the presence of X-ray emission 
from this remarkable object. The brightest region of X-rays 
from this source is clearly extended, seems to sit within the 
shell of emission defined by the "head" of the radio source, and 
has a highly absorbed spectrum consistent with being either an 
extremely steep power law or a collisionally ionized plasma at 
a temperature ^0.6 keV. 

We have considered the possibility of the Tornado being 
powered by jets from an X-ray binary, of corresponding to the 
synchrotron nebula powered by an energetic rotation-powered 
pulsar, or of representing an unusually-shaped supernova rem- 
nant. While the available data do not allow us to definitively 
rule out any of these alternatives, the spectrum and morphology 
of the X-ray emission seen by Chandra, when combined with 
evidence at other wavebands for the interaction of a shock with 
a molecular cloud, all suggest that the head of the Tornado is 
a "mixed morphology" supernova remnant, producing thermal 
X-rays from the interior of a limb-brightened radio shell, and 
expanding into a dense and possibly complicated distribution 
of ambient gas. 

This argument provides a reasonably simple interpretation 
for the "head" of the Tornado, but the bizarre morphology of 
the "tail" remains unexplained. We speculate that the Tornado 
might result from expansion of a SNR into a highly elongated 
progenitor wind bubble, but explanations involving pre-existing 
helical magnetic fields, or high-speed outflows from a still un- 
seen central source, can also not be discounted. Deep X-ray 
observations with Chandra and XMM are needed to better de- 
termine the morphology and spectrum of X-rays from the head, 
and to study the possible X-ray emission from the tail. 

We thank Dan Harris, Jasmina Lazendic and Fred Seward for 
useful discussions, and Crystal Brogan for supplying us with 
her radio images of the Tornado. This work was supported by 
NASA through SAO grant GO2-3041X (BMG) and by the NSF 
(JMM). 
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Table 1 



Spectral fits to X-ray emission from the head of the Tornado. 



Model 


N H (10 zz cm"- 1 ) 


T / kT (keV) 


/ x (l(r u ergcm- z s" 1 ) 


xi/y 


RS 




6 +0 - 8 

U '°-0.3 


-5.1 


10.2/9= 1.1 


PL 


12* 


>4.5 


> 1.2 


10.8/9= 1.2 


PL 


-3.4 


2.3 (fixed) a 


-0.05 


20.3/10 = 2.0 


PL 


5±1 


3.5 (fixed) a 


-0.2 


15.7/10= 1.6 



"See §§3.2 and 3.3 for details. 

Uncertainties and lower limits are all at 90% confidence. All models assume interstellar absorption using the cross-sections of Balucihska-Church & McCammon 
(1992), assuming solar abundances. Models used: "PL" indicates a power law of the form f e oc e~ r where T is the photon index; "RS" indicates a Raymond-Smith 
spectrum of temperature T (Raymond & Smith 1977). Fluxes are for the energy range 0.5-10 keV, and have been corrected for interstellar absorption. 
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17 h 41 m s 50 s 40 s 30 s 20 s 10 s 



RA (J2000) 

FIG. 1. — Radio and X-ray images of G357.7-0.1 ("the Tornado"). The top panel shows a 1.4-GHz VLA image of the Tornado, corrected for primary beam 
attenuation and with a spatial resolution of 14" X 11" (Brogan & Goss 2003). (The "eye" of the Tornado, is an unrelated Hll region; Brogan & Goss 2003.) The 
image in the bottom panel shows Chandra ACIS-I data in the energy range 0.3-10 keV after exposure correction and adaptive smoothing. In the bottom panel, black 
contours represent X-ray emission at levels of 76%, 84% and 92% of the peak brightness of region A, while white contours represent radio emission at levels of 7.5 
and 50 mJy beam -1 . Point sources identified with wavelet ect but which are too faint to see in the smoothed image are marked with "+" symbols. The position of 
the OH maser is marked with a "X" symbol. X-ray emission seen in the northeast and southeast corners of the image is due to the large exposure correction near the 
edges of the field of view, and does not correspond to real emission. Regions discussed in the text are indicated. 



